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Abstract 

We propose a new method that fixes the coupling to two photons of 
the recently found lightest QCD resonance, the a. This coupling provides 
crucial information for discriminating the yet unknown nature of this spe- 
cial state. Our method uses available data on the nucleon polarizabilities 
together with analyticity and unitarity. Taking into account all the uncer- 
tainties, our result is r po i e = 1.2 ± 0.4 keV. 
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1 Introduction 



The lowest isospin 1 = and angular momentum J = QCD resonance is usually 
called the a and plays a special role in the QCD dynamics and in the QCD non- 
perturbative vacuum structure. Recently, it has been shown that the a is also the 
lowest QCD resonance by fixing the mass and width of this state with a precision 
of just tens of MeV in Ref. [1]. Making an analytic continuation of the 1 = 
and J = partial wave S-matrix in the region of validity of Roy equations, these 
authors find a zero at E — [(441+g 6 - % (272l? 2 )] MeV on the first Riemann sheet 
which reflects the a pole on the second Riemann sheet. Also on the first Riemann 
sheet, the inverse of the partial wave irir S-matrix S — 1 + 2i f3(t) T{t) has a zero 
at E*. Here, 

Tit) = W) cot[5(t)] +/?(*)' (1) 

5(t) is the scalar-isoscalar tctt phase shift, (5{t) = ^Jl — Am 2 / 1 and t = E 2 . The 
position of the a resonance pole has been confirmed in [2] at E — [(484 ± 17 — 
% (255 ± 10)] MeV. 

The nature of the a remains one of the most intriguing and difficult issues 
in particle physics. There are have been many proposals about its substructure: 
q — q state, it — tt molecule, (qq) — (qq) tetraquark, glueball, and of course, several 
admixtures of these substructures. The size of a — > 77 width can shed light on 
this question. 

77 — > (7T7r)/ =0 2 amplitudes have been calculated using twice-subtracted dis- 
persion relations in [3] and [4] including the recent data on tx-k final state inter- 
actions which contain the a pole in the scalar-isoscalar contribution. They get 
4.09 ± 0.29 keV [3] and 1.68 ± 0.15 keV [4] for the a into two photons width. 
The origin of this discrepancy is discussed in [4]. More recently, the authors of 
[5] made an amplitude analysis of the world published data on 77 — > 7r + 7r~ and 
find two regions of solutions. The width o — > 77 in these regions are predicted 
in Ref. [5] to be 3.1 ± 0.5 and 2.4 ± 0.4 keV, respectively. 

2 Method 

In Ref. [6] , we proposed a new method that fixes the coupling to two photons g ail 
of the a meson found in the nn scattering amplitude [1, 2] using only available 
precise experimental data on the nucleon electromagnetic polarizabilities together 
with analyticity and unitarity. This differs from the analysis of [7], where the 
properties of the a meson of a Nambu-Jona-Lasinio model were used. Nucleon 
electric a and magnetic (3 polarizabilities are well measured using low energy 
Compton scattering on protons and neutrons with a + /3 constrained by a forward 
dispersion relation [8]. The results are [9]: a cxpt = 12.0 0.6 and j3 cxpt = 1.9 =F 0.5 
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for protons and ct expt = 11.6 1.5 and (3 expt = 3.7 =F 2.0 for neutrons. Here and in 
the rest of the paper, polarizabilities are given in 10" 4 fm 3 units. 

The authors of Ref. [10, 11] wrote a sum rule for a — (3 using a backward 
dispersion relation for the physical spin-averaged amplitude. The s-channel part 
of this sum rule is related to the multipole content of the total photo-absorption 
cross-section. While the t-channel part is related through a dispersion relation 
to the imaginary part of the amplitude which using unitarity is given by the 
processes 77 — > tttt and nn — > NN [11]. The result is the BEFT sum rule, 

1 f°° du I u~ 

{a -0 = 2? L ^f +2 if r ^ = yes) - ° (A * = ,io)| 

+ £«fri^{l/Mira*)! 

where M p is the proton mass, the partial wave helicity amplitudes /+(£) and 
f+(t) for NN — > mi are Frazer and Fulco's [12], and the partial wave helicity 
amplitudes F^t) and F^t) for 77 — > tttt were defined in [13]. The s-channel 
part of the integrand is obtained from that of the forward physical amplitude 
by changing the sign of the non-parity flip multipoles (An = no) and yields 
(a - p) s = -(5.0 ± 1.0) [14]. The "experimental" (a - is therefore 15.1 ± 1.3 
for protons and 12.9 ± 2.7 for neutrons. Products of helicity amplitudes in Eq. 
(2) appear only as moduli products, which can be negative if the phases differ 
from the mi phase shift in an odd number of 7r's. The d-wave contribution to 
(a — PY is much smaller than the s-wave one; hence, it is a good approximation 
to take just the Born result (a — ~ —1.7 [15]. Finally, the "experimental" 
value for (a — /3)q is 16.8 ± 1.3 for protons and 14.6 ± 2.7 for neutrons. The input 
\Fq (t)\ for (a — (3)q in Eq. (2) is what we want to fix using this "experimental" 
value. The other input for this quantity, the Frazer- Fulco's [/+(£) I amplitude, 
is known with enough accuracy for our purposes from the old determination in 
[16], though it could be improved using recent data on irn phase shift. We have 
assigned a 20 % to the theoretical determination of (a — /3)q from this source. 
Notice that the 1/t 2 factor in the integrand of (a — /3)q makes the well known 
low-energy, and to a lesser extent, intermediate-energy contributions to be the 
dominant ones. 

On the physical sheet, we can write the twice-subtracted dispersion relation 
[17], 

f °(o = Lit) - m 

where c is a subtraction constant fixed by chiral perturbation theory (CHPT) 
[17, 18], c = a/AMfl with a ~ 1/137 the fine-structure constant, f w = 92 A MeV 



ct^ 



e_ y°° dt' Ljt'ymVL- 1 ^') 
7i J Ami t' 2 t' — t — ie 



(3) 
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the pion decay constant, 



Q(t) = exp 



dt' 5(f) 



7T J Ami t' V — t — IS 



(4) 



is the scalar-isoscalar tttt Omnes function [19] which gives the correct right-hand 
cut contribution and L(t) is the left-hand cut contribution. In this way we ensure 
unitarity, the correct analytic structure of Fq (t) and that the a pole properties 
enter through the scalar-isoscalar phase-shift 8{t) from T(t) in (1). Here, we shall 
use a simple analytic expression for T(t), compatible with Roy's equations, which 
takes a three parameter fit from [2] including both low energy kaon data and high 
energy data. This fit is valid up to values of t of the order of 1 GeV 2 , which is 
enough for the (a — (3)q integrand in Eq. (2). 

At the a pole position on the first Riemann sheet [3, 4] 

where e is the electron charge, g% n7r is the residue of the nn scattering amplitude 
at the a pole on the second Riemann sheet and g all g a7T7T is proportional to the 
residue of the 77 — > scalar-isoscalar scattering amplitude on the second Rie- 
mann sheet. The proportionality factors are such that g an7T and g ail agree with 
those used in [3]. The pole width is given by 

r ( \ "l^J (p\ 

r po ie(^ -> 77) = — 4M ^ (6) 

that agrees, modulo normalizations, with that of Ref . [4] . 

Low's theorem fixes the amplitude F^t) to be the Born term at very low- 
energy [20]. As first approximation, we therefore consider the left-hand cut con- 
tribution L(t) in (3) to be the Born contribution L B (t), 

This leads to Fq^Ib which when input in the sum rule (2) gives (a — 
P)o\b — 6.7 ± 1.2, where the uncertainty comes mainly from the input data 
on I . This result is 5a away the "experimental" values quoted above. 

When analytically continued to complex t, the amplitude Fq(£)| b has a pole at 
t a = {[(474 ± 6) - i(254 ± 4)] MeV} 2 with g a7TW = [(452 ± 4) + i(224 ± 2)] MeV 
and using (5) we get gtm/g^B = (0.49±g:gf) - i (0.37 ± 0.03) which leads to 
r po ie(o' — ► 77) |s = 2.5 ± 0.2 keV. But this Fo(i)|s does not reproduce the "ex- 
perimental" (a — PY and hence we need to go beyond the Born approximation 
L B (t) for the left-hand cut. 
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Figure 1: The integrand of (a— (3)$ in (2). The dashed line is when using L(t) = Lsit) 
in (3) and the continuous line is when using L(t) = Lsif) + L A (t) + Ly(t) in (3) as 
explained in the text. 



The first corrections to Ls{t) originate in the resonance exchange 77? — > R — > 
77r, with R = cti, p and u [4, 17]. In the narrow width approximation, the a\ 
exchange contribution to L(t) is 



L A {t) 



C 



327T/2 



M 2 

t+ m los 



' l+/3(t)+t A /t * 
1 - 0(t) + t A jt j 



(8) 



while the p and u resonances exchange contribution to L(t) in nonet symmetry 
(M p = M UJ = M V ~ 782 MeV) is 



L v {t) 



3 v 



t 



My_ 
Pit) 



log 



l-P(t)+t v /t / 



(9) 



with tR = 2(M|, — m 2 ). The low energy limit of Ly(t) goes as t 2 and we fix R v = 
1.49 GeV -2 from the well known uj — > tcj decay. Though the low energy limit of 
L A (t) goes as t and corresponds to the charged pion electromagnetic polarizability 
(a — /3) w ± or equivalently to L 9 + L w = (1.4±0.3) • 1CT 3 in CHPT [21], we consider 
L A (t) as an effective contribution for moderate higher values of t with C a real 
constant to be determined phenomenologically and not connected to the pion 
polarizability. This is supported by the fact that the ain^ interaction is not so 
well known at intermediate energies. We fix the effective C by requiring that the 
"experimental" value of (a — /3)q is reproduced within 1.5 standard deviations of 
the total uncertainty when L(t) in (3) is given by L(t) — Lb{£) + L A (t) + Ly{t). 
This procedure leads to C — 0.59 ± 0.20 and the integrand of the sum rule is 
given in Fig. 1 as a continuous line. Notice that the zero at to in the integrand 
of (a — P) Q in (2) when using F® (t)\ B has clearly disappeared now. 

The low-energy 77 — > 7r°7r° cross-sections obtained when the left-hand cut 
is either Ls(t) or the full L(t) case studied above are similar [4]. The central 
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values are compatible with data for values of t below (450 MeV) 2 and are above 
data but compatible within two standard deviations for larger values of t up to 
(600 MeV) 2 and within one standard deviation for t between (600 MeV) 2 and 
(800 MeV) 2 . 

3 Results and Conclusions 

The scalar-isoscalar amplitude F^t), obtained using L(t) = L B (t)+LA(t)+L v (t) 
as explained above, is analytically continued to the complex plane, and at t a 
on the first Riemann sheet one gets g ayy /g a7T7T = (0.23l ;og) — « (0.30 ± 0.03) 
which has a smaller absolute value when compared with g cll j ga^B and leads to 
r p oie(°" —> 77) = (1-0 ± 0.3) keV. The error quoted here is from the uncertainties 
in the "experimental" value of (a — (3)q and the inputs of the sum rule (2) only. To 
obtain the rest of the uncertainty, we vary the a properties in the 7nr scattering 
as follows. We still use the three parameter fit formula including low energy 
kaon data and high energy data for cot(8(t)) in [2] as input in the amplitude 
T(t) but with parameter values slightly modified to reproduce the a pole position 
t a = ([(441±6)-i (272±4)] MeV) 2 found in [1]. In this case, we get g a7T7r = [(480± 
7) + % (191 ±3)] MeV. With that T(t) input in the dressed Born amplitude in Eq. 
(3), one gets (a - (3)l\ B = 6.1 ± 1.1, g cll l 9<t^\b = (0.57 ± 0.02) - i(0.41 ± 0.03) 
and r pole (o" — > 77) | £ = (3.8 ± 0.4) keV. The integrand of (a — (3)q in (2) for this 
case is very similar to the dashed line of Fig. 1. The effective value of C in (8) 
is C = 0.62 ± 0.20 when fixed to reproduce the "experimental" value of (a — (3)q 
within 1.5 standard deviations of the total uncertainty. With this new C, the 
analytic continuation to t a gives g ail j 1 g a ^ — (0-3lio!o?) — * (0.32 ± 0.03) and 
r po ie(o' — > 77) = (1-5 ± 0.4) keV. Again, the integrand of (a — (3)q in (2) for this 
case is very similar to the continuous curve of Fig. 1. 

Making a weighted average of the two cases discussed above, we get 

r po ic(^ -> 77) = (1-2 ± 0.4) keV . (10) 
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